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ABSTRACT 

We use a collection of 14 well-measured neutron star masses to strengthen the case that a substantial fraction 
of these neutron stars was formed via electron-capture supernovae (SNe) as opposed to Fe-core collapse SNe. 
The e-capture SNe are characterized by lower resultant gravitational masses and smaller natal kicks, leading 
to lower orbital eccentricities when the e-capture SN has led to the formation of the second neutron star in a 
binary system. Based on the measured masses and eccentricities, we identify four neutron stars, which have 
a mean post-collapse gravitational mass of ^1.25 M , as the product of e-capture SNe. We associate the 
remaining ten neutron stars, which have a mean mass of ~4.35 M Q , with Fe-core collapse SNe. If the e- 
capture supernova occurs during the formation of the first neutron star, then this should substantially increase 
the formation probability for double neutron stars, given that more systems will remain bound with the smaller 
kicks. However, this does not appear to be the case for any of the observed systems, and we discuss possible 
reasons for this. 

Subject headings: stars: neutron, stars: evolution 



1. INTRODUCTION 

Precise neutron star mass determinations, coupled with a 
theoretical knowledge of the pre-collaps e mass, can be used 
to test the neutron star equation of state. Podsiadlowski et al. 
(2005) used the double pulsar system J0737-3039 for such a 



test. They inferred from the low mass (1.249 ± 0.001 M Q ) 
of Pulsar B in this system, that it formed via an electron- 
capture supernova ( SN), an inference first made in|Podsiad- 
lowski ~et al.| ([2004) and independently by |van den HeuveT 
( |2004[ ). In addition to such uses, well-measured neutron star 
masses are extremely helpful in understanding the formation 
scenarios of these objects. 

In this work, we consider the entire sample of known neu- 
tron stars with well-measured masses. We find that the mass 
distribution is most compatible with the existence of two dis- 
tinct populations, a higher-mass (^1.35M Q ) and a lower- 
mass (^1.25M Q ) population. We interpret these two popu- 
lations to be the result of distinct evolutionary formation sce- 
narios: the low-mass population originates in electron-capture 
SNe and has received low kicks, while the high-mass popula- 
tion is the result of iron core collapse SNe. 

In §2 we compare and contrast the two principal channels 
for the production of neutron stars: e-capture supernovae and 
Fe core-collapse supernovae. The current sample of 14 well- 
measured neutron stars is presented and discussed in §3; these 
all have mass uncertainties of < 0.025 M . In §4 we perform 
some statistical tests which provide support for the hypothe- 
sis that there are two parent populations of pre-collapse core 
masses. We carry out a simple population synthesis study 
in §5 of the expected eccentricity distributions for e-capture 
and Fe core-collapse SNe, using their different anticipated 
core masses and natal kick speed distributions. In §6, we at- 
tempt to fit the observed systems with well determined neu- 
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tron star masses into the two principal evolutionary scenarios: 
the "standard" channel and the double-core channel. We sum- 
marize our results and draw some general conclusions in §7. 
In particular, we find that (i) a substantial fraction of neutron 
stars are formed in e-capture SNe, and (ii) there is evidence 
from our work that the double-core formation scenario is less 
unlikely than previously thought by most workers in the field. 

2. EVOLUTIONARY HISTORY 

Neutron stars are believed to form through two main evolu- 
tionary channels: iron core collapse and electron-capture su- 
pernovae. The first occurs in a massive star when it has devel- 
oped an iron core which exceeds the Chandrasekhar mass and 
no more nuclear burning can take place. The resulting mass of 
the neutron star depends not only on the neutron star equation 
of state, but also on the mass of the iron core and the maxi- 
mum iron core mass for which a successful supernova can oc- 
cur. The latter depends on the details of the supernova mech- 
anism that are still not fully understood. If the iron core mass 
is too large, the explosion mechanism fails and the core col- 
lapses to a black hole. In the presently most popul ar paradigm 
of delayed neutrino-driven ex plosions (see, e.g.,|Mezzacappa 
et al.|2007 Janka et al. 2008} , the explosion takes place when 
enough neutrino energy has been deposited in the gain region 
outside the proto-neutron star to overcome the binding energy 
of the remaining core, stop the accretion and initiate an out- 
flow. The characteristic energy of such a delayed explosion 
has to be of the order of the characteristic binding energy of 
the remaining core (~ 10 51 ergs). Hence, successful iron-core 
collapse supernovae are expected to have explosion energies 
close to this characteristic energy. 

In contrast, an electron-capture (e-capture) supernova oc- 
curs in a very degenerate ONeMg core, long before an iron 
core has developed, and is triggered by the sudden capture of 
electrons onto Ne nuclei, taking away the hydrostatic support 
provided by the degenerate electrons (e.g., Nomoto 1984). 
This occurs at a cha racteristic density (~ 4.5 x 10 y gcm" J ; 
Podsiadlo wski et al.|2005| l, which in turn can be related to a 
critical pre-collapse mass for the ONeMg core of ~ 1.37M . 
Hence, an e-capture supernova is expected to occur when 
a degenerate ONeMg core reaches this critical mass either 
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Table 1 

Comparison of Fe-Core Collapse and e-Capture Supernovae 



Properties 



Iron Core Collapse 



e-Capture SN 



Supernova Properties 
Explosion energy 
Ejecta 



Neutron Star Properties 
Masses 

Neutron star kick 



Binary Properties 
Occurrence 

Eccentricity 
Recycled pulsar spin 



~ 10 51 ergs 

rich in heavy elements 

(Fe, Si, O) 



range of masses 

large standard kick 
(cr ~ 265 km s-') b 



single or binaries 
high 

mis-aligned with orbit 
(e.g. geodetic precession) 



< 10 50 ergs a 

few heavy elements 



characteristic mass 
~ 1.25M 
low kick 



preferentially 
in binaries c 
low 

aligned with orbit 



Dessart et al. 1 2006 1; Kitaura et al. 1 2006 1 
Hobbs et al.H27305T 
Podsiadlowski et al.|{2004) 



by accretion from an envelope (in side an AGB star [e.g., 
Siess||2007| |Poelarends et al.| [2008 1 or in a helium star [e.g., 
Nomoto||1987[ ), by accretion from a c ompanion star (so- 
called accretion-induced collapse [e.g., Nomoto & Kondo 
1991 1), or as a consequence of the merger of two CO white 
dwarf s and the su bsequent formation of an ONeMg core (e.g., 
|Nomoto & Iben|1985[ l. Since the collapse occurs at a charac- 
teristic ONeMg core mass, the resulting neutron-star mass is 
entirely determined by the equation of state and the amount 
of core material that is ejected in the supemovan The case of 
Pulsar B in the double pulsar system J0703-3039suggests that 
this mass is close to 1 .25 M Q (Podsiadlowski et al. 2005). Fur- 
thermore, since essentially the whole core collapses to form a 
neutron star, the remaining envelope is relatively easy to eject, 
leading to a fainter super nova with the ejection of very few 
heavy elements (see, e.g., |Dessart et al.||2006{ |Kitaura et al.| 
2006 ). It has recently been argued that the large kicks most 
neutron stars receive at birth (Hobbs et al. 2005) are caused 
by an accretion shock instability that causes a wobbling of 



the core, imparting momentum in the process (e.g., Blondin & 



Mezza cappa 2006] |2007[ |Foglizzo et al.|2007 but see [Fryer 
& Young ( 2007 ) for a more skeptical point or view). Since, 



<5i Young ( 2007 )Tor a more skeptical point < 
in the case of an e-capture supernova, the explosion occurs 
before these instabilities have time to grow, no large kick is 
expected for a neutron star formed through this channel. 

The suggestion that electron-capture supernovae may pro- 
duce low supernova kicks and a distinct low-mass n eutron 
star populatio n was firs t made independently b y 
lowski et al.| ( |2004|> and |van den Heuvel| (|2OO4ji0 



'odsiad- 



van den 



Heuvel| ( |2004| i specifically discussed this low-mass, low-kick 
population in the context of binary radio pulsars, and used 
the then-current observations of several NS-NS binaries and a 
NS-WD binary to argue that they formed via e-capture. 

Table 1 summarizes the main differences in neutron-star 
and supernova properties for these two channels. Note, in 

4 This ignores the role of rotation which may be important, in particular, 
in the case of an accretion- or merger-induced collapse. 

5 The latter author also suggested a third more massive population of neu- 
tron stars with masses around 1.85 Mq from stars with an initial mass around 
20 M Q . 



particular, that, for neutron stars formed from iron-core col- 
lapse, one expects a range of masses that is determined by 
the range of iron core masses in the progenitors that allows a 
successful explosion, while in the case of neutron stars from 
an e-capture supernova one expects a fairly well determined 
mass. Thus, the distribution of post-supemova neutron star 
masses directly constrains not only the equation of state, but 
also the properties of successful supernova explosions. 

3. NEUTRON STAR SAMPLE 

There are 14 neutron stars which have masses known with 
an accuracy of better than ~ 0.025 M . The majority of these 
(twelve) are from double neutron star systems; two are in bi- 
nary systems with suspected white dwarf companions. The 
properties of thes e systems are summarized in Table[2](for ref- 
erences see, e.g., |Stairs|2008[ ). A histogram of the measured 
gravitational masses is shown in the top panel of Figure [T] 
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Figure 1. Mass histograms for the sample of 14 neutron stars. (Top Panel): 
The measured (gravitational) masses of the neutron stars. (Middle Panel): 
The masses of the neutron stars corrected for accretion as discussed in the 
text. (Bottom Panel): The precollapse (baryonic) masses of the neutron stars, 
based on one particular illustrative neutron-star equation of state. 



The rapidly rotating pulsars have likely been spun up by 
the accretion of a small to modest amount of matter (0.001 - 
O.O7M ). We correct for this effect by subtracting the mass 
which would be necessary to spin up the star, treating it as a 
classical uniform-density sphere accreting from a disk that ex- 
tends down to its surface. We have verified that for a range of 
plausible equations of state for neutron star matter, more so- 
phisticated treatments lead to accreted (gravitational) masses 
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Table 2 

14 Well Measured Neutron Star Masses 



Pulsar Name 


Mass of Recycled 


Mass of Young 




Eccentricity 


Pulse Period 


Reference 




Neutron Star (Mq) 


Neutron Star (Mq) 


(hours) 




(ms) 




J0737-3039A/B 


1.3381 ±0.0007 


1.2489 ±0.0007 


2.4 


0.088 


23 


[Kramer et al. (2006J 


B1534+12 


1.3332±0.0010 


1.3452 ±0.0010 


10.1 


0.273 


38 


|Stairsetal.||2TO2y 


J1756-2251 


1.32 ±0.02 


1.24 ±0.02 


7.67 


0.18 


28 


JsTarrs|j20l)8| 


J 1906+0746 


1.365 ±0.018 


1.248 ±0.018 


3.98 


0.085 


144t 


|KasIan|pD08} 


B1913+16 


1.4414 ±0.0002 


1.3867 ±0.0002 


7.92 


0.617 


59 


|Weisberg & Taylorjl|2005} 


B2127+11C 


1.358±0.010 


1.354±0.010 


8.05 


0.681 


30 


Ijacobv et al. 120051 


J 1909-3744 


1.438 ±0.024 


white dwarf 


36.7 


< io~ 6 


2.9 


Jacoby et al. 12005) 


Jl 141-6545 


white dwarf 


1.27±0.01 


4.74 


0.172 


393 f 


|Bhat efaI7|poD8| 



Note. — All known neutron stars with a mass measured with better than 0.025 Mq accuracy. 
t These periods are said to be associated with the "young pulsar". 



that differ from our simple model by less than 



10% (see, 

The results are shown in the mid- 



[T] Note the high degree of similarity of 
tha 



e.g., Cook et al. 
die panel of Figure 

this histogram with tKat for the uncorrected masses; the max- 
imum mass correction for any one neutron star is ^0.07 M Q 
(for J1909-3744). The corrections for the other neutron stars 
were less than ^0.02 Mq. 

Finally, we used a repre sentative equation of state for neu- 
tron star matter ("MPA", Miither, Praka sh, & Ainsworth|1987| l 
to translate the observed gravitational mass into a pre-collapse 
mass by calculating the baryonic mass corresponding to each 
gravitational mass. The results are shown in the bottom panel 
of Figure [T] In general, the pre-collapse masses are shifted 
upward by ^0.13 Mq. 

As the equation of state remains theoretically uncertain, we 
calculated t he corrections for each of t he equations of state 
collected in Lattimer & Prakash (2001). Within this collec- 
tion, the correction to the mass of a 1.25 Mq neutron star var- 
ied over the range 0.09-0.18 Mq. However, given the small 
range in mass considered (1.25 - 1.4 Mq), the choice of equa- 
tion of state has little effect on the relative correction between 
any two systems within this range. The net result of choos- 
ing a different equation of state would be a systematic shift 
in the bottom panel of Figure [T| as opposed to any significant 
stretching or skewing. 

One can see from Figure [T]that there are two apparent pop- 
ulations of neutron star mass: one centered at ~1.25Mq and 
~1.35M© (post collapse) and ~1.37M Q and ~1.48Af (pre- 
collapse - for the assumed "MPA" neutron-star equation of 
state). The higher of these two mass groups is suggestive of 
an origin in an Fe core collapse SN, while the lower of the 
two groups likely comes from electron-capture SN events. 

4. STATISTICAL TESTS 

We make use of two statistical tests to try to quantitatively 
evaluate our hypotheses, the Kologorov-S mirnov (KS) tes t 
and the Anderson-Darling (AD) test (e.g., |Press et al.|20Q7| >. 
The AD test is more powerful as it takes into account the in- 
tegrated difference between the cumulative distributions one 
is comparing, while the KS test considers only the maximum 
difference. 

The first test we perform is for normality, checking whether 
the distribution is consistent with a single Gaussian which has 
the mean and standard deviation of the observed populations. 
The cumulative distribution for the observed neutron stars is 
shown in Figure [2] as filled circles connected by a dashed his- 
togram. The cumulative distribution for the single Gaussian 
described by a mean mass of 1 .325 M Q and standard deviation 
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Figure 2. The cumulative distribution of neutron star masses. The observed 
distribution is shown by the black dots and dashed black line. The CDF for 
the best two population model is shown by the solid red line (see text). The 
blue curve is the CDF for a single Gaussian with the population mean and 
standard deviation. 



of 0.056M Q is shown as the blue curve in Figure [2] We were 
not able to reject this hypothesis of a single Gaussian mass 
distribution with a KS test, but were able to marginally reject 
it at the 70% confidence level with the AD test. 

We also tested the hypothesis that there are two distribu- 
tions present, each of which is represented by a Gaussian. The 
best such fit is given by: 



dN/dM = 0.707 e 



(M-1.345) 2 /2cr, 



+ 0.293 e 



(M-1.246)-/2(T2 



(1) 



with al = 0.025 ± 0.004M Q and a 2 = 0.008 ± 0.005 Mq. The 
uncertainties on the mean masses of the two Gaussians are 
1.345±0.003M o and 1.246 ± 0.003 Mq, respectively. The 
cumulative distribution function corresponding to this distri- 
bution is shown plotted as a solid red curve in Figure [2] The 
amplitudes of the two Gaussians are, as expected, reflective 
of the fact that 4 of the 14 neutron stars are in the lower-mass 
group. As can be seen from the cumulative distribution for 
the double Gaussian, the fit is very good, compared to a single 
Gaussian with the population mean and standard deviation. 

While these statistical tests and fits do not, by themselves, 
constitute a proof of two populations, coupled with the other 
pieces of evidence (i.e., appropriate system eccentricities and 
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Figure 3. Comparison of system eccentricities of the double neutron star bi- 
naries. (Top Panel) : The measured masses of the young neutron star plotted 
against the system eccentricity. We identify the low-mass, low-eccentricity 
systems as being the result of e-capture SNe, marking them as red. An 
amount 0.01 was artificially added to the eccentricity of J0737 to separate 
it from J1906 on the plot. (Bottom Panel): The results of the Monte Carlo ec- 
centricity simulation described in the text. The black curve is the distribution 
of Fe-core collapse systems and the red curve is for electron-capture systems. 

theoretically expected masses for e-capture SNe), they do lend 
support for the hypothesis of two populations. 

5. ECCENTRICITY CALCULATIONS 

In order to illustrate what the eccentricity distributions of 
e-capture vs. Fe core collapse SNe might look like, we car- 
ried out the following simple statistical study. In all cases, we 
assume that the second SN explosion takes place with a He 
or CO core in a circular orbit with the first-born neutron star. 
We take the core mass to be in the range 1.5-2.OM for the 
e-capture scena rio and 2.5-6.0Mq for the Fe core collaps e 
scenario (e.g., Podsiadlowski et al. 2004} [Dewi et al.| |2005 ). 
The orbit is assumed to have been circularized during a prior 
episode when the evolving core expands sufficiently to trans- 
fer at least a small amount of mass to the first-born neutron 
star, thereby spinning it up to millisecond rotation periods. 
The orbital separation at the time of the second SN explo- 
sion is taken to be uniformly distributed over the range 1-3 
times the orbital separation needed for the He or CO core to 
fill its Roche lobe (a result expected from more detailed pop- 
ulation synthesis calculations, e.g., Dewi et al. 20061. Finally, 
the natal kick distribution is taken to be a Maxwellian with 
a equal to either 30 km s" 1 or 265 km s" 1 , for th e e-capture 
and Fe core collapse scenarios, respectively (e.g., De wi et aT] 
|2006[ l. The resulting orbit al eccentricities are computed from 



the expression given by Brandt & Podsiadlowski ( 1995). The 
above assumptions should hold for both the standard forma- 
tion channel as well as for the double-core channel. 

For each scenario, the initial system parameters (i.e., core 
mass and orbital separation) and the natal kick were chosen 
via Monte Carlo techniques for 10 7 systems. The results are 
shown in Figure [3] The red curve represents the eccentricity 
distribution for the e-capture scenario while the black curve 
is for the Fe core collapse explosions. Note that in the lat- 
ter case the eccentricity distribution is rather broad (in fact 
some ~60% of the systems become unbound) while for the 
e-capture events the eccentricity distribution peaks at e ~ 0.2, 
and extends down to rather low values of e. In fact, we ob- 
serve that these two distinct distributions are rather consistent, 
respectively, with the three cases we tentatively identify with 
e-capture SNe and those three with higher eccentricities that 



we associate with Fe core collapse SNe. 

6. CONSISTENCY OF THE SCENARIOS 

There are two main scenarios for the formation of double 
neutron star systems. They both start with a pair of massive 
primordial stars, i.e., with masses between ^8 and 25 M Q . 
In the "standard scenario" the more massive star evolves first, 
fills its Roche lobe, and stable, quasi-conservative mass trans- 
fer to the secondary may occur if the mass ratio is not too 
extreme and if the initial orbital period is in the right range 
(i.e., ~ a month to a year). In this scenario, the system can 
not undergo a common envelope phase during the first stage 
of mass transfer, otherwise the orbit would not be sufficiently 
wide after the formation of the first neutron star to allow for 
the successful production of the second neutron star (as dis- 
cussed below). Given the fact that the orbit should initially 
be wide, it is actually advantageous for the first supernova ex- 
plosion to occur via an e-capture SN, in order to yield a small 
natal kick and thereby help the system remain bound. There 
is indeed a class of such wide binaries containing a neutron 
star in nearly circular orbit with a massive d onor (see, e.g., 
|Pfahl et al.|[2002l |Podsiadlowski et al.||2004] . After the for- 
mation of the first neutron star, the companion star evolves, 
fills its Roche lobe, and the subsequent mass transfer onto 
the neutron star leads to a common envelope phase due to the 
extreme mass ratio of the system. In order for the common 
envelope phase to avoid a merger of the neutron star and the 
core of its companion, the orbit must be wide to begin with, 
as alluded to above. The result of a 'successful' common en- 
velope phase is a He or CO core in close (< 1 day) orbit with 
the first neutron star. The original neutron star is spun up via 
accretion from its He/CO star companion as it evolves toward 
core collapse. If the subsequent supernova explosio n is via 
e-cap ture (for He-star companions of mass < 2M Q ;[Nomoto] 
1984), then the natal kick may be small and the final binary 
pair of neutron stars would have a modestly sm all eccentric - 
ity and low systemic space velocity (|van den Heuvel|[2007 1. 
Whether the second neutron star would form via an e-capture 
SN or Fe core-collapse SN would depend on the original mass 
of the secondary star and the orbital period after the first neu- 
tron star has been formed. 

Both this standard scenario and the "double core" scenario 
that we describe next are summarized in schematic diagrams 



in |Podsiadlowski et al.| ( 2005] Fig. 1) and |Dewi etaL| ( [2006 
Fig. 1) The various possible orderings for the different core- 
collapse mechanisms for the two neutron stars in the standard 
scenario are summarized in Table J3] along with some quali- 
tative comments on the likelihoodof each. In particular, in 
a recent detailed population synthesis study of DNSs, which 
included both the standard and the double-core channel, Bel- 



czynski et al. ( 2009) predicted that category II systems should 
be by far the most dominant systems. From the above discus- 
sion and the notes in Table [3] we would expect double neutron 
stars formed via the standard scenario to consist of a recycled 
pulsar with a mass indicative of an e-capture SN formation 
(Categories II and IV; i.e., ~ 1.25M ). Moreover, if the or- 
bital eccentricities are low, then both neutron stars would most 
likely be formed via e-capture SNe (i.e., Category IV). Nei- 
ther of these expectations is borne out by the observational 

data in Table 2. 

In the "double-core" scena rio (Brown 1995' |Bethe &] 
Brown|[T998 Dewi et aL||2006| >, the primordial binary is re- 



quired to have a pair of stars whose mass is the same to within 
3 - 7% of each other. Given that massive stars seem to of- 
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Table 3 

Order of Fe-Core Collapse vs. e-Capture Supernovae 



Category Neutron Star Formation Type and Order Standard Scenario Double Core Scenario Observed? 



I 


Fe core collapse + Fe core collapse 


possible 


probable 


yes 


II 


e-capture + Fe core collapse 


most favored 


inconsistent 


no 


III 


Fe core collapse + e-capture 


possible 


probable 


yes 


IV 


e-capture + e-capture 


possible 


some fine tuning 


no 



ten reside in binaries whose stars have comparable mass, this 
~5% "window" is not nearly as rare as it might seem. In fact, 
roughly speaking some several percent of all massive stars 
may occur in such comparable-mass binaries if the binary 
fraction is high and if the mass ratio distribution is roughly 
flat. If the two stars have comparable mass, and the orbit is 
relatively wide (i.e., ^months to years) then both stars will en- 
ter a double common-envelope phase once the primary starts 
to transfer mass to the secondary. At that point the primary 
is expected to have evolved a CO core while the secondary 
will also be evolved, but will more likely have only a He core. 
Both of these cores will spiral in inside the common envelope 
formed of the envelopes of both stars. The result will be a 
close pair of a CO and a He core (or less likely a pair of CO 
cores). The orbital period of the primordial binary should be 
wide (i.e., months to years) if the cores are to avoid merger 
during the common-envelope phase. The CO core evolves 
first, most likely, to an Fe core-collapse supernova and leaves 
a ~1.35 Mq neutron star (i.e., Categories I and III; see Ta- 
ble]^]), though it could also experience an electron-capture SN 
(Category IV). When the He core evolves, it expands some- 
what and transfers some mass to the first neutron star, thereby 
recycling it. If the He core mass is relatively low (< 2M ; 
Nomoto 1984), it will go on to undergo an e-capture super- 
nova, leaving a lower-mass neutron star with a smaller natal 
kick (i.e., Category III). Thus, in this scenario, the recycled 
pulsar is the more massive of the two, having been formed in 
an Fe core-collapse supernova, while the second, an unrecy- 
cled pulsar, should be lower in mass. 

This double-core scenario is consistent with the data for the 
six neutron star binaries. Three of them have higher orbital 
eccentricities, e = 0.27,0.62,0.68, with both neutron stars in 
the system having masses consistent with Fe core-collapse su- 
pernovae (Category I). The other three have lower eccentric- 
ities, e = 0.08,0.09,0.18), and the unrecycled pulsar is con- 
sistent with having undergone an e-capture supe rnova (with 
mass ~1.25M ), i.e., Category III (see also van den Heuvel 



2007). The eccentricities and NS masses of these six systems 
are plotted in Figure [3] Since we have only six systems, it is 
premature to draw any firm conclusions, though it is remark- 
able that none of the observed systems falls into the a priori 
most favored standard scenario category. 

In the case of the two neutron star-white dwarf binaries, the 
scenario is likely somewhat different than described above. 
One of the neutron star-white dwarf systems (J 1141-6545) 
has a significant eccentricity of 0. 17 and a lower-mass neutron 
star. This suggests that, in this case, the white dwarf actually 
formed before the neutron star, otherwise the second phase 
of mass transfer, leading to the formation of the white dwarf, 
would have circularized the orbit. This is the likely conse- 
quence of a first phase of conservative mass transfer where the 



initial masses of the binary components were relatively close 
and the primary had a mass just below the minimum main- 
sequence mass for neutron-star formation (~ 7M ). Af- 
ter this first mass-transfer phase, the secondary has accreted 
enough mass to end its evolution as a neutron star (see, e.g., 
|Church eT"a l. 2006). In addition, after the ensuing common 
envelope phase, the secondary is more likely to have a mass 
just above the minimum neutron star formation mass and this 
naturally favours an e-capture collapse. The formation of this 
system via e-captur e was previously suggested by |van den| 
|Heuvel| ( |2004| |2006). The other neutron star-white dwarf sys- 
tem (J 1909-3744) is highly circularized and it is likely that the 
neutron star formed first and the orbit was circularized dur- 
ing a common envelope phase involving the progenitor of the 
white dwarf. 

7. SUMMARY AND CONCLUSIONS 

We have shown that the population of 14 well-measured 
neutron-star masses is consistent with being comprised of a 
subset of four that were likely the result of electron-capture 
SNe, while the others resulted from Fe-core collapse SNe. 
The lower neutron star masses (^1.25M ) of the four candi- 
date electron-capture SNe and relatively low orbital eccentric- 
ities (<? < 0.18) of the systems that contain them are both good 
indicators of this type of formation mechanism. The remain- 
ing 10 neutron stars have larger masses (^1.35M ) more in- 
dicative of Fe core-collapse SNe; in the systems where these 
more massive neutron stars formed second, larger orbital ec- 
centricies are observed, consistent with a larger natal kick. 

We discussed four categories of formation scenarios for 
producing neutron stars in close binaries, especially double 
neutron stars. These include the possibilities that (i) either the 
first or second neutron star was formed via an Fe core collapse 
SN and (ii) that either neutron star might have been formed in 
an electron-capture SN. Any of these four possibilities could 
be connected with either the "standard scenario", in which a 
common envelope phase occurs only after the first SN explo- 
sion, or with the so-called "double-core" scenario, wherein 
both comparably massed primordial stars are simultaneously 
stripped of their envelopes. (See Table|3]for a summary of the 
eight possible combinations.) 

None of the observed systems falls into Categories II or IV, 
because the recycled pulsar is never the lower-mass product of 
an e-capture supernovae. For producing double neutron star 
systems, it appears that the standard scenario is somewhat dis- 
favored. The evidence that we have presented and examined 
seems to favor Category I (both the young and the recycled 
NS are higher-mass) and Category III (a recycled higher-mass 
NS and a young lower-mass NS), both forming via the double- 
core channel (see Table [3]). 

If our interpretation of the neutron star mass data is correct, 
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i.e., that the double-core channel is the preferred scenario for 
producing double neutron stars, this has profound implica- 
tions for the fate of neutron stars within a common envelope. 
Recall that the original motivation for proposing the double- 
core channel was the estimation, by |Chevalierj ( 1993| |1996[ l 
that a neutron star within a common envelope might undergo 
hypercritical accretion and collapse to a black hole. On the 
other hand, there is evidence that at least some neutron stars 
survived a common envelope without b eing conve rted into a 
black hole (e.g., PSR B0655+64, |Tauris et al.|2000) . And, as 
noted above, the neutron star in j 1909-3744 also appears to 
have spiraled through a massive common envelope. However, 
it seems possible that systems with massive white-dwarfs sec- 
ondaries could also be produced via the double-core channel. 
We note that, even though the double-core channel requires 
very special initial conditions, such systems are more compact 
when the first supernova occurs and are therefore much more 
likely to survive as bound systems after the first supernova 
than in the standard channel. As a consequence, the birthrate 
of double neutron star systems in the double-core channel can 
be comparable to the stan dard model: using binary popula- 
tion synthesis simulations, |Dewi et al.| (2006 ) estimated their 
birthrate to be 1 0~ 6 - 10~ 5 yr" 1 (bu t with substantial uncertain- 
ties, see Belczynski et al. ( |2006| l for another estimate). This 
is probably sufficient to account for the observed number of 



double neutron-star systems (see, e.g., Kalogera et al.|2004) . 
Nonetheless, whether or not neutron stars can survive a com- 
mon envelope, we believe we have provided some further sup- 
port for the double-core channel formation of double neutron 
star systems. 

We thank Chris Fryer, Michael Kramer, Onno Pols, and Ed 
van den Heuvel for extremely helpful discussions. 
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